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Abstract 


It is shown that an indirect synthesis method can be used in the 
efficient optimal design of multi-degree of freedom, multi-design 
element, nonlinear, transient systems. The technique begins with a 
limiting performance analysis which requires linear programming for a 
kinematically linear system, following which the system is selected 
using system identification methods such that the designed system re- 
sponds as closely as possible to the limiting performance. The ef- 
ficiency is a result of the method avoiding the repetitive systems 
analyses accompanying other numerical optimization methods. 
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1. Introduction 

An essential ingredient of the conventional methods in the optimum 
design of dynamic systems has been an iterative technique that progres- 
sively searches for the open design parameters which optimize some index 
of system performance without violating the imposed constraints. While 
this search procedure presents no serious problems for simple systems 
with few design parameters, its practicality diminishes rapidly with 
increasing system complexity, and becomes formidable in terms of compu- 
tational requirements for a large, real world system. The main difficulty 
therein stems from the fact that for each iteration, the entire system 
dynamics must be solved and that convergence to a global minimum, if it 
can be reached, usually requires a large number of iterations. 

Recently a new approach called the "indirect synthesis method" 
has been suggested (Ref. 1) and successfully demonstrated for a single 
degree of freedom system. Satisfactory results have also been obtained 
by applying this approach to a nonlinear impact absorber (Ref, 2) and 
to vehicle suspension systems (Ref. 3). The basis of the new approach 
is to select the open design parameters such that the response of the 
portion of the system being designed approximates the limiting performance 
solution which is independent of the hardware devices used. It suffices to 
solve the system dynamics only once. As a consequence, the new approach 
offers great computational advantages over the conventional methods. 

Also, by carrying out the limiting performance solution for several 
constraint levels it yields the limiting performance characteristics 



(trade— off curve) of the system under consideration. Apart from its 
use in automated design, this information, which is not available from 
conventional sources, is of considerable value to the mechanical designer 
in providing the feasibility characteristics of his proposed design^ 

Details of the formulation of the new design approach are presented 
in Ref 1, Essentially the method consists of two distinct phases: 
time optimal synthesis and parameter identification. One first replaces 
the design elements, regardless of configuration, by generic or control 
forces u(t). The limiting performance solution of the resulting system 
is sought in the time domain of interest, yielding the optimal (ideal) 
generic forces u*(t) and other ideal system responses. In the second 
phase of the design one selects, element by element, the design para- 
meters that respond most closely (in an optimal sense) to the limiting 
performance solution. The elements being designed can be linear, non- 
linear, passive or active. The limiting performance problem is particularly 
simple when those portions of the system not being designed and the over- 
all dynamics are linear, for then the problem can be formulated as one 
of linear programming. The element by element parameter design can be 
achieved by unconstrained curve fitting techniques such as least square 
or min-max fits. Thus the nonlinear programming optimization problem 
normally encountered in the direct methods is now replaced by a more 
efficient linear programming problem. 

As stated earlier, the new technique has been applied to single 
degree of freedom systems only. Also, previously the control forces - 
were represented in the limiting performance problem by piecewise 
constant bursts of forces (of a preselected time interval) over the time 
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period of concern. This often results in a response that reverses its 
direction suddenly and repeatedly, thus making it difficult to 
identify parameters by curve fitting. In order to ease the task of 
parameter identification, we have reformulated the transient limiting 
performance problem using "smooth’* control forces, i.e, ones repre~ 
sented by Fourier series. It is the purpose of this work to outline the 
formulation and to demonstrate the new methodology for multi-degree, 
multi-design element systems subject to transient disturbances . 

2, Limiting Performance Formulation 

Consider the problem of designing portions of a dynamic system 
subject to transient disturbances so that some index of performance is 
minimized (or maximized) and certain response constraints are satisfied. 
Replace these portions, which can be nonlinear in the state variables, 
by control forces u(t) and suppose that the remaining portions of the 
system are linear as are the overall kinematics. It is also assumed 
that both the performance index and constraints are linear functions of 
the state variables. For such systems, the limiting performance problem 
can be formulated as one of linear programming. 

2 . 1 Equations of Motion 

The general equations of motion for the system described above 


can be written in the form 



where 


X = displacement vector 

u = control force vector which replaces portions of the system 
being designed 
f = forcing function vector 

M = N X N mass matrix, N being the number of degrees of freedom (DOF) 
C = N X N damping matrix 
K = N X N stiffness matrix 

_V = N X J coefficient matrix associated with control force 
vector, J being the number of control forces 
= N x L coefficient matrix associated with forcing function 
vector, L being the number of forcing functions 


We choose to represent the control forces in terms of "smooth” 
functions, e.g. Fourier series. Previously (Ref, 1) "nonsmooth" piece- 
wise constant control forces are employed. It will be demonstrated that 
the smooth control forces lead to better designed systems. Each control 
force will be represented by a finite-term Fourier series of preselected 
basic frequency: 


NF 


u = y (a. cosw t + b, sxnu) t) + a.^, j - 1>2,. , . , J 
i in n jn n ju 

n=l 


m 

i 

n=l 


( 2 ) 


= y (A cosw t + B sinu) t) + A, 
^ ^ n n n n i 


or u 
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where 



and o) = nit/T, T being the half period. Note that this representation 
n 

of the control forces places no linearity restrictions on them. 

Similarly, the relative displacements are represented by 
NF 

X. = y (p, cosu) t + q. sino) t) + p.rk» ^ ” 1*2,. . , N 
i ^ - ^in n ^in n iO 

n=l 

_ NF _ 

or X = y (P cos(o t + Q sino) t) + P. (3) 

n n n n 0 

n=l 

where 



and, 

_ NF __ 

f - y (R cosio t + S sinoj t) + R_ 
^ n n n n 0 

n 


( 4 ) 
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with 


R = 
n 


In 


2n 


Ln 


S = 
n 


In 


2n 


Ln 


,and Rq = 


10 


20 


LO 


Differentiation of (3) results in 


= y (- 0 ) P srnoD t + 0 ) Q cosoj t) 
^ nn n n n n 
n 


(5) 


and X = y (- 0 )^ P cosw t - Q sinw t) 
n n n n n n 
n 


C6) 


Now substitute equations (2) - (6) into (1) . This yields 


y M (-03^ P COSO) t - Q sinu t) + y C (-u) P sina> t + w Q cosw t) 
-nn n nn n ^—nn n nn n 
n n 


+ y K (P cosu) t + Q sino) t) + K P^ = y V (A cosw t + B sinw t) 
^— n n n n — 0 n n n n 

n 


+ V A_ + y F (R COSO) t + S sinw t)+ F R 


n 


Collecting terms of cosoi^t, sinoj^t and constants on both sides and making 
use of the orthogonality on (-T,T), we obtain 


(-0)2 M + K)P + 0) CQ = VA + F R 
n — — n n— n — n — n 


-0) C P + (-o)^M + K)Q = VB + F S 
ir-n nr— — n — n — n 


( 7 ) 


KP„ = VAq + F Rq 


for n = 1,2,... ,NF. 
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Eqns. (7) are recognized to be a set of simultaneous linear equations 
in P^, and P^, the right hand side being the Fourier coefficients of 

u and the forcing functions. Thus we have shown through Eq. (7) 
that the Fourier coefficients of x are related to those of u. 

Frequently, e.g., in the case of many shock isolation systems, 
the equations of motion can be written as 

^ = J/u + Ff (8) 

For this special case, it is necessary that we express the relative 
displacement x in terms of the Fourier coefficients of the u’s by direct 
integration of (8) instead of solving Eqn. (7). Thus, one integration 
of (8) gives 

NF 

^ = la 

n 


sinoint _ cosojnt _ - _ 

“ ^ B^) + vA^t + dt +Ci, 


n 


(On 


(9) 


and after a second integration 


^ COS(0j,t _ 

Mx = y V( - 

- n ^ ^ 


sinojnt 

— 2 S 

n 


) + + +C,t + C2, (10) 


0 2 


where Ci and C 2 are integration constants to be determined by the 
initial conditions, 

2.2 Limiting Performance 

With the state variables now expressed linearly in terms of the 

Fourier coefficients A-, A , B of the control forces u(t), and f being 

U n n 

a known quantity, it is a simple matter to express the performance index 
and the constraints linearly also in terms of these coefficients. 
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The limiting performance problem theti becomes: 

Find the Fourier coefficients A^, A , B , where n=l,2,. . .,NF, 

0 n n ’ ’ 

such that the performance index 

ip = max max|(j)^(AQ, A^, B^, f, t) | , i=l,2,, , ,,I 
is minimized subject to the constraints 

* k=l,2,* . , ,K 

where the are I response functions and are the upper and 

lower bounds of the constraint 6^. This problem is readily recognized 

as one of linear programming (LP) and it can easily be reduced to the 
standard LP form (Ref . 1) . Then any available LP codes can be used 
to determine these coefficients and hence the optimal responses u*(t), 
X*(t), x*(t), x*(t). 

3. Results 

The new design approach with "smooth" control force representation 
has been applied to a number of transient systems. Two selected 
examples are presented below. 

3.1 TWo Degree of Freedom System 

Consider the two degree of freedom system shown in Fig. 1(a). 

A step velocity of Vq “ 1 in/sec is applied to the rigid base over a 
period of 1,5 sec. The design problem is to determine the kj^, k 2 > 
and C2 so that the peak acceleration of is minimized subject to the 
constraints that the rattlespaces and X 2 be less than 3/4 in. 




(b) 


^^8’ ^ System with isolators replaced by 

control forces. 
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The substitute system for the limiting performance study is shown in 
Fig. 1(b). The governing equations are 

V 

Zi + ui = 0 
Z2 “ ni + U2 = 0 

with the kinematic conditions 
X 2 = - y 

The limiting performance problem then is to minimize 
= max] I 

with the constraints 

1x^1 <3/4 
IX 2 I 1 3/4 

The limiting performance solutions obtained by a linear programming 
technique are shown in Fig, 2 and 3 together with the results of para- 
meter identification obtained by fitting u*(t) and u*(t) respectively to 
+ c^x*(t) and k^x^Ct) + c^x^Ct) . 

Here for simplicity, an ordinary unconstrained least square 
technique is used. The curves labeled "parameter identification" are 
the solutions to the equations of motion of the designed system. The 
limiting performance solution yields a peak acceleration of 1.35 for 
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m]^ while the designed system gives 1,60. It is also seen that the 
rattlespaces are slightly violated, by 8 % for the upper one and 5.5% 
for the lower. 

For purposes of comparison, the same problem was solved using the 
piecewise constant representation of the control forces (Ref . 1) * 

The results are shown in Fig, 4 and 5. It is apparent that although 
the limiting performance characteristic of maxj zj | — max] ] are almost 
identical for the two representations, the task of curve fitting the 
isolator responses in this case (Fig. 4(b) and 5(b)) is much more 
difficult than the case of "smooth” representation. Here the worst 
rattlespace (upper) violation is as much as 100 % and the performance 
index is also considerably higher than the "smooth" method. 

3.2 Six Parameter Designs of Three Degree of Freedom System 

In order to make quantitative comparisons of the new design ap- 
proach with conventional computer optimization, the structure of Fig, 6 (a) 
was studied. The goal is to find the parameters k 2 , k 3 , C 2 , 03 

such that for base inputs y = t^e ^ the maximum of the peak values of 
I ^2 1 » 1^1 is minimized subject to 

]x^(t) + aG(t)|j< £ 1 , |xi(t) - e0(t)l <_ £i 
|x2(t) -b6(t)l G2» |x2(t) +d0(t)l ±^2 




Fig. 4 Upper and lower isolator response histories - force-time 
(by piecewise constant representation) . 



LOWER RATTLESPACE Xj UPPER RATTLESPACE 
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Fig. 5 Upper and lower isolator response histories - displacement 
time (by piecewise constant representation) . 
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Fig, 6 (a) Three degree of freedom system. Initial conditions: 

z^(0) - z^(0) = 0(0) = z^(0) = z^(0) = 0(0) = 0. 

(b) Substitute system used for limiting performance study. 
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The substitute structure employed in the design process is shown in 
Fig. 6(b). The equations of motion are 

z = -(ui + U4) 

Z 2 " -('^2 + U 3 ) 

e- = -—(ui + U4) + ^ u - ^ U3, 

with the kinematic relations 

xi = zi - yi, X2 = Z2 - Y2 

and the initial conditions 

zjCO) = Z 2 ( 0 ) = zi(0) = ^ 

The radius of gyration is designated by p. 

The design parameters are determined by matching the ideal isolator 
forces uj(t), i=l,2,3, to 

gl = - be*) + -^ ( xj *- be* ) 

g 2 = ^ ( xj + de*) + -^ ( x5 + de* ) 

g 3 = ki(xi + a9*) + (x* + a9* ), 

For a=b~d=25, p=23,5> ei=e2=0*l>a least square curve 

fitting technique leads to = 2.994, k 2 = 7.129, k 3 = 1.626, cj = 2.716, 
C 2 = 5.062 and C 3 = 0.426. The active objective functions are plotted in 
Fig. 7 together with the ideal isolator forces. It is seen that the 
peak objective function is 0.83 for the designed system, compared to the 
value of 0.6 given by the limiting performance solution. The constraint 
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Fig. 8 Constraint responses for the designed system. 
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responses are plotted in. Fig. 8. The worst violation of the constraints 
is about 25% . For the limiting performance study we chose NF“2 in the 
Fourier expansion and a half period of two seconds was selected. This 
is close to the time at which the base velocity drops to zero and its 
displacement reaches its maximum. 

This particular design problem was chosen because the same system 
was designed using classical techniques in Ref. A. There the publicly 
available (ASIAC, Ref. 5) optimization program POSI was employed. This 
standard computational optimization program selected a system with an 
optimal objective function of 1.15 versus our 0.83. Our more optimal 
design, of course, is partly a result of constraint violations. This 
apparent problem, which is considered below, was suggested in Ref. 3 as 
being so serious as to make questionable the value of this method. 

Perhaps the most interesting result is the comparison of computer 
times. Typical designs (such as the one just presented) made by the 
new approach took about one to two percent of the computer processor 
time expended by the classical computer optimization using POSI. 

The problem of constraint violation is not as serious as it may 
seem. Techniques more sophisticated than unconstrained least square 
fits, e.g., a constrained min-max fit, may be used to reduce the constraint 
violations. Constrained parameter identification techniques are also 
useful if it is known that for practical consideration the design para- 
meters must be restricted to certain values. In addition, it is possi- 
ble to correct constraint violations by finding the design parameters 
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based on matching the response being violated (Ref. 1). Also, this 
problem can be remedied by using our optimal design parameters as 
starting values in a conventional computer optimization program. 

For the above problem a design that satisfied all constraints was 
obtained after one iteration by POSI. 

4. Conclusions 

It has been demonstrated that the indirect synthesis design 
approach can be used to obtain a near optimal design for multi-degree 
of freedom and multi— design— element systems - at least for the inputs 
considered here. The design occasionally violates the desired con- 
straint levels. As explained, however, corrective measures can be 


taken . 
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